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Abstract. An extended Hubbard model with phonons is considered on a D- dimensional 
lattice. The symmetries of the model are studied in various cases. It is shown that for a certain 
choice of the parameters a superconducting SU q (2) holds as a true quantum symmetry — but 
only for D = 1. 

1. Introduction 

In this article we are going to study extensions of the Hubbard model on a D- 
dimcnsional lattice and their symmetries. 

The Hubbard model was originally introduced in [1] and is the simplest model de- 
scribing systems of itinerant interacting electrons in solid-state physics. Its importance is 
mainly due to the possibility that it may describe high-temperature superconductivity. 

Since the work of Yang and Zhang [2,3] it has been known that the Hubbard model 
has a SU(2) x SU (2))/Z2-symmetry. This symmetry is the product of two separate SU (2) 
symmetries: a magnetic symmetry, which accounts for the (antiferro-)magnetic properties 
of the electron system, and a superconductive symmetry, which is supposed to give rise 
to superconductivity when it is broken. 

A. Montorsi and M. Rasetti investigated whether the symmetry of the standard 
Hubbard model generalizes to a quantum group symmetry of the Hubbard model with 
phonons. In [4] they claim the existence of a superconductive S'C/ 9 (2)-quantum group sym- 
metry and a standard magnetic SU (2)-symmetry for a Hubbard model with non-local 
phonon interaction. We were able to verify this symmetry for a particular extension of 
the Hubbard model, but find it to be restricted to 1-dimensional lattices. Moreover, this 
Quantum Symmetric Hubbard Model is not the standard Hubbard Model with Phonons. 
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2. The Extended Hubbard Model 

The Hamiltonian of the extended Hubbard model is given by [4,5]: 

TT tj(Ioc) . tt(Ioc) tt(Ioc) . Tj(non-loc) f-l\ 

where 

H ei° C) = U X! n ^ ni i ~ A* X! nilJ ^ 
i i,a 

i v ' 

^-l = -A-E(^T+«aK (4) 

ff (no»-M = ^^ T .. a t ffa . CT (5) 

with T y - = Tt given by 

Ty = te- h ^- R exp{Cfly • - f ,•) + i£ • - ft)} (6) 

In these expressions the Hamiltonian depends on the operators , a ia , ft , x*i , where 
djcr, aj CT are fermionic annihilation/creation operators for an electron with spin a = {j, J.} 
at site z, while ft, ir, are momentum and displacement operators for an ion at site i with 
the usual commutation relations. 

The Hamiltonian also depends on R\ which gives the rest position of the the ion at 

site i. Note that the expression Ru = j-^ — iri always has unit modulus and that in the 

J | Hi — iri j | 

one-dimensional case it just amounts to a sign. \R\ — Rj\ is the interatomic distance at 
equilibrium so that it does not depend on i, j. In the following we will explain the various 
terms composing Hnub'- 

H^° c ^ is the local electron-electron interaction and is the sum of two terms. The first 
one, the on-site repulsion u^2 i ni^riii, is reminiscent of the Coulomb repulsion between 
the electrons and is determined by the parameter u. The second one, the chemical poten- 
tial —A* Si a n * CT > shows that the Hamiltonian is written in a grand-canonical formalism; 
the parameter fi fixes the avarage number of electrons in the lattice. 

Hph^ i s the kinetic term for the phonons, which are described by a set of decoupled 
Einstein oscillators, all with the same frequency u and the same mass M. 

H-ti-yh is the local phonon-electron interaction term [6,7] with coupling constant A, 
which at each site gives an attractive force proportional to the number of electrons and 
to the ion displacement. 

The crucial term is the non-local one. It is a hopping amplitude, and gives the prob- 
ability that an electron can jump from one site to another one. Notice that we have 
retained only the nearest neighbour terms (ij) and hence assumed negligible overlap be- 
tween all other atomic orbitals. In the extended model this amplitude depends on the ion 
displacement and momentum, resulting in a non-local interaction between the phonons 
and the ions. 

Let' s study different special cases of the extended Hubbard model and see which 
models can be recovered. 
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2.1 Overview of the different models contained in the extended Hubbard model. 

In the following table we give a list of different models contained in the Extended 
Hubbard Model. A more detailed explanation follows in the remarks. 

I) Hubbard Model with: ( = 

a) generic case: A ^ 0, k ^ 

. . t* . Lang-Firsov-transform. . 

b) with local phonon interaction: A ^ 0, k — < ==>■ c) 

c) no local phonon interaction: A = 0, k ^ Mca n field ggp rox - standard Hubbard model 

d) Standard Hubbard Model plus decoupled Einstein Oscillators: A = 0, k = 
II) Hubbard model with: £ ^ 0, but k = 

a) with local phonon interaction: A ^ 

b) no local phonon interaction: A = 
III) General case: C^0, /? ^ 0, A ^ 

Remarks 

ad I) When k = 0, ( = 0, but A ^ (case lb) the Hamiltonian Huub can be used to describe 
bipolarons [8], a model in which there is only the local electron-phonon interaction. 
Notice that equivalently, it is possible to describe the bipolarons with an Hamiltonian 
of the type Huub with A = , but k/0 (case Ic) . This can be done by performing a 
unitary transformation, the Lang-Firsov transformation [9], on , a \ a , pi , Xi with a 
unitary operator 

U(k) = exp ^in ■ Y^Pi n i<^j ■ ( 7 ) 

Performing a mean-field approximation on the phonon variables (when A = 0), one 
recovers the standard Hubbard model, 
ad II) When k = the Hamiltonian Hnub describes the Hubbard model with phonons 
added. To see this, remember that the hopping amplitude in (5) is originally defined 
by Hubbard [1] as 



Th = J d D r t>*{r-R t -Xi) {^^j ^(f-Rj-Xj). 



(8) 



where ^(r — Ri) is the Wannier electron wave function. is a function only of 
dij = (Ri + Xi) — (Rj + Xj). We approximate the Wannier electron functions with 
atomic orbitals, which show an asymptotic exponential decay ^>(r) ~ e - *^ and find 



*(f). (9) 



Because of the rapid exponential decay of \&(r), we can neglect r in \r — o^l so that 
[10] 



V 3ii T(3y) = -C^-T(a l3 ) (10) 



I a ij I 



which integrates to T(a tj ) = T e fl a «l. 
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\aij | = \Ri — Rj + Xi — Xj | can be approximated by \£i — x 3 \ <C \Ri~ Rj | such that [5] 

T.^te^f-C ^-y ^-xj)) (11) 
\ \Ri-Rj\ ) 

with a new constant t = T cxp(—(\Ri — Rj\). 

3. Symmetries of the Extended Hubbard Model 

None of the terms added to Hubbard model affect the magnetic S'?7(2)-symmetry. 
Hence, in the sequel we shall only be concerned with the superconductive symmetry of 
the model. Denote the Jimbo-DrinfePd generators of U q (su(2)) with X ± ,H. They satisfy 
the following commutation relations (see e.g. [11]): 

-,H „-H 



\H,X& 



= ±2X< ± \ \X+,X-]=1 q -r- (12) 

q-q 1 



3.1 Local commutation relations 

Definition Local representation of the superconductive U q (su(2)) at each site I [4]: 

p s {X+)=K^ =e- i **a\,a\ i (13) 

p s (X-) ee <"> = e^'a nai1 = (tf <+>)t (14) 

p s {H)=2K\ z) =n IT +n u -l (15) 

The parameter 4> appearing in (13), (14) does not affect the [7 g (su(2))-commutation 
relations, and for the moment it should be regarded as a free variable, which will be 
determined by the commutation relations with the Hamiltonian. ^ can be interpreted as 
the parameter of a Lang-Firsov transformation [9] of the fermionic operators. 

Theorem 1. The local part of the Hamiltonian commutes with the local generators 
if and only if the following conditions are satisfied: 

H = - - -Mu 2 h 2 <f> 2 = — — — — (18) 
' 2 4 2 Muj 2 y ' 

3.2 Global commutation relations 

Switching signs on p s (X ± ) gives again a representation of U q (su(2j). It is necessary 
to consider both representations, i.e. 

pf(X+) = ±p s (X+), pf(X-) = ±p s (X-), pf(H) = Ps (H), (19) 

if the convention is chosen that fermionic operators at different sites anticommute. For 
each lattice site I a sign a(l) G {1,-1} and the associated representation pl^ will be 
determined by the symmetry. The local commutation relations are not affected by this 
choice. 

Further, it is necessary to fix some ordering of the lattice sites to be able to define a 
tensor product and hence to construct a global symmetry. 
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Definition Global representation of the superconductive U q (su(2)) [4]: 

k {+) = (g) P f i) (a^- 1 ) (x+))=j2 <^(o n eQ ^' ^ (+) n e_Q *^ > ' ( 2 °) 

l l r<l r>l 

K { - ] =<g)p^(M N - 1 \X-)) = 5>(0]I eQ *^' II e~ aK r\ (21) 

z z r<Z r>l 

JfM =0p^)(A( JV - 1 )(F)) = ^^ W , (22) 
z z 

A is the coproduct of U q (su{2)) (see e.g. [11]), a = ln(g) is the deformation parameter, 
is the number of lattice sites. 

Theorem 2. The non-local part of the Hamiltonian commutes with the global gener- 
ators 

J^(+) jj(non-loc) 

if and only if the following conditions are satisfied for i,j nearest neighbours: 



a(i) = -a(j), (24) 

2/? = $, (25) 

Rea = -Rij ■ $(h, (26) 

J] e aK ^ } = [] r- a * K ^- (27) 

i<r<j i<r<j 



Remarks to Theorem 2: Eq. (24) imposes that nearest neighbours must have 
opposite signs. This gives a restriction on the possible lattices, e.g. a triangular lattice is 
inconsistent with our conventions. 

Eq. (25) fixes the parameter $ of (13), (14) and relates it to the parameter k appearing 
in the Hamiltonian Hu-ub (6). As k can be interpreted as the parameter of a Lang- 
Firsov transformation on the Hamiltonian, while ^ can be interpreted as the parameter 
of a Lang-Firsov transformation on the generators of the symmetry, eq. (25) can be 
interpreted as a consistency relation, requiring the same transformation to be done on 
the Hamiltonian and on the generators. 

Eq. (26) determines the real part of the deformation parameter a of the quantum 
group. In order for such equation to make sense it is necessary to choose a lexicographic 
ordering of the lattice sites, so that the sign of the term does not depend on the 
particular couple i, j. Notice that there is no restriction on the imaginary part of the 
deformation parameter; we can safely choose it to be real. 

Eq. (27) is trivially satisfied if a = 0; it is a strong condition in the case a ^ 0: It then 
implies for nearest neighbours i < j that there cannot be a site r in between them, i.e. 
with i < r < j. However, such a condition requires that the lattice on which the Hubbard 
model is defined is one-dimensional, and that a ordering of sites is chosen in which the 
sites are numbered from left to right in increasing or decreasing order. 

Proof. (Sketch) The (undeformed) generator commutes with fj( non - loc ) given 
by (5) as in the classical (q = 1) case. 



k { -\h^ 



-loc) 



(23) 
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We still have to calculate 



-HR ir R) (J. J. J. J 



( a li a h a ]i a n a ^ a Ji ~ a ii a li a ]i a ]i a n a j-[) 



(i<j) 



Zij [ 2 cosli ( ^Rij ■ $(hj -2 cosh Q#ij • + 



( 2 cosh ( ^Rij ■ $(ft^ - 2cosh Q-Ry • $C?i+ \^ 



Z y -(e5 a * - 1) + Zji ^ e -Aj-*Cft e -ia _ ^ 
Zy(e*y*<M Q ' - 1) + Z^e"^ - 1) 



x 



x: e ^)c-^-»e-*ftatot n ea ^' n e ~ Q ^ 

( {i,j),i<Kj r<l,r=£i r>l,r=j£j 



where we have introduced the abbreviation 



n e 

r>i,r^j 



(28) 



The sums containing different numbers of fcrmionic operators are linearly independent 
and must all vanish separately. (We have fixed a normal ordering for the fermionic opera- 
tors.) Studying the term containing aj^atj — a^at^, the conditions (24)-(27) are obtained. 
It turns out that the same conditions guarantee that also the other sums vanish. 

3.3 Summary 



Classical Symmetry 

a = 

D arbitrary 


Quantum Symmetry 

a ^ 0, q = cxp(a) 
D = 1 


Id), lib): S?7(2) s -symmetry; fj, = § 

la): S'C/(2) s -symmetry; 

H = u/2- h 2 K 2 Muj 2 , A = hMuj 2 K 

Ib),Ic),IIa): no symmetry 


III): S'/7 g (2) s -symmetry; 
Rea = 2k(U, 

li = u/2- Ti 2 k 2 Muj 2 , A = hMuj 2 K 
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Remark: After choosing an appropriate ordering of the lattice sites for model III) 
in the 1-dimcnsional case, we have set k = —Rij • it. The models Ib),Ic) and Ha) have no 
symmetry, because the condition 2k = j-j^jt is not satisfied. 

4. Discussion 

The standard Hubbard model Id) (with decoupled Einstein oscillators) and the Hub- 
bard model lib) without local phonon interaction (A = 0) have a superconductive SU(2) S 
symmetry at "half-filling" /i = u/2. Two other models la), III) with local phonon inter- 
action (A 7^ 0) have superconducting symmetries SU(2) S and SU q (2) s respectively at 
(i = u/2 — h 2 K 2 Muj 2 . It is well known [9] that a Lang-Firsov transformation can be 
performed on the model la), shifting the parameters and thereby eliminating the local 
phonon interaction A. The resulting Hamiltonian has a S[/(2) s -symmetry at "half-filling" 
/i = u/2. It turns out that a Lang-Firsov transformation with the unitary operator (7) can 
also be used to eliminate the local electron-phonon interaction term from the Extended 
Hubbard Model. The parameters are shifted according to 

A -> A - Muj 2 M, u -> u - 2h\ ■ k + Mu?h 2 K 2 , fj, -> \i + h\ ■ k - \/2Mu?n 2 R 2 (29) 
and the hopping term H^°^ h loc ^ becomes 

t e^i^-^qh a) a a ia (l + (g"4 - l)m,-*) (l + («* - lK_ ff ) + h.c. (30) 

(i<j>,<7 

The resulting Hamiltonian has a superconducting SU q (2) s quantum symmetry at "half 
filling" n = u/2 (!) with A = and only for D = 1. 

It is even possible to eliminate all phonon terms from this Quantum Symmetric Hamil- 
tonian by a mean field approximation [8,9] without breaking the quantum symmetry. More 
precisely one performs a thermal average over i/p/j-eigenstates and assumes uncorrelated 
Einstein oscillators (note that A = !). The exponential in the hopping term is then 
approximated by a temperature-dependent constant. 
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und Umwelt for financial support. 
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